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Simple scalar field cosmological models are considered describing gravity assisted crossing of the 
l/") \ phantom divide line. This crossing or (de)-phantomization characterized by the change of the sign 

of the kinetic term of the scalar field is smooth and driven dynamically by the Einstein equations. 
Different cosmological scenarios, including the phantom phase of matter are sketched. 

(N ; 

^| PACS numbers: 98.80.Cq, 98.80.Jk 

. I. INTRODUCTION 

m ; 

The discovery of cosmic acceleration [jj has become one of the main challenges for modern theoretical physics, 
^v^j \ because it has put forward the problem of the so called dark energy, responsible for this acceleration. This dark 
. energy should possess a negative pressure such that the relation between pressure and energy density is less than 
f-^. ' —1/3 0- A considerable amount of dark energy models was proposed (see, e.g. which now have to confront 

OO , with growing quantity of observational data. 

Being enigmatic by itself, the phenomenon of the cosmic acceleration - dark energy - negative pressure hides 
inside another subphenomenon, which is not yet so well confirmed experimentally, but if confirmed, it would become, 
perhaps, more intriguing than all the preceeding cosmological puzzles. We mean here, the phenomenon of the so called 
phantom dark energy, i.e. the matter with the relation between the pressure and dark energy k = p/ p less than —10. 
Some observations indicate that the present day value of the parameter k < — 1 provides the best fit. According to 
some authors, the analysis of observations, permits to specify the existence of the moment when the universe changes 
the value of the parameter k from that the region fc > — 1 to fc < — 1 @. This transition is called "the crossing of 
tp.' the phantom divide line". Other authors speak about double crossing of the phantom divide line (jj]])- In any case 
the existence of the phantom phase of cosmological evolution put forward some important theoretical problems. For 
example, considering a perfect fluid with barotropic equation of state, where k = const and k < — 1 one can easily see 
that a universe during its evolution after a finite amount of time encounters a special kind of cosmological singularity, 
characterized by an infinite value of the velocity of expansion, which has been dubbed "Big Rip" 0, [ij • Being a little 
bit unusual in comparison with well accepted Big Bang and Big Crunch singularities, the Big Rip singularity is quite 
similar to them from the formal (geometrical) point of view. Notice , th at in the modern cosmological literature also 
other "exotic" types of singularities are under consideration [nj, HH UH UM • 

Another difficulty connected with the problem of phantom energy has a more fundamental character. It is well- 
known that the models with minimally coupled scalar field are very efficient both for the description of the inflationary 
stage of the development of the early universe and for the explanation of the late-time cosmic acceleration phenomenon 
However, it is easy to see, that the standard minimally coupled scalar field cannot give rise to the phantom dark 
energy, because in this model the absolute value of energy density is always greater than that of pressure, i.e. |fc| < 1. 
A possible way out of this situation is the consideration of the scalar field models with the negative kinetic term. 
Again, being a little unusual from the point of view of common wisdom, these models have become quite popular 
novadays. 

Perhaps, the most thrilling problem arising in connection with the phantom energy is the crossing of the phantom 
divide line. The general belief is that this crossing is not admissible in simple minimally coupled models and expla- 
nation of this phenomenon requieres more complicated models such a multifield ones or models with non-minimal 
coupling between scalar field and gravity (see e.g. pH llH ITU llS^. 

In the present paper, we would like to explore the consequences of the possibility of existence of simple minimally 
coupled scalar field models allowing the crossing of the phantom divide line. As a matter of fact, we would like to 
present some toy examples, suggesting that such a crossing could be unavoidable and dictated by the continuity of the 
solution of Einstein equations. In the next section we discuss in some detail a simple toy model, giving an illustration 
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of the change of sign of the kinetic term in the scalar field cosmology. In the third section we sketch some other 
models and in the last section we attempt to discuss the above consideration in a more general context. 

II. TOY MODEL WITH A CROSSING OF THE PHANTOM DIVIDE LINE 

Let us consider a simple cosmological model, representing a flat Friedmann universe with a metric 

ds 2 = dt 2 ~ a 2 (t)dl 2 (1) 
filled with a barotropic fluid with an equation of state 

P = hp, (2) 

where p and p are energy density and pressure of the fluid respectively and — 1 < k < 1. Choosing conveniently the 
normalisation of the constants in the theory one can write the Friedmann equation in a simple form 

h 2 = p, (3) 

where the Hubble variable is defined as h = - . The energy conservation law is 

p=-3h(p + p) (4) 

which together with Eq. J2J gives immediately 
From Eqs. J2J and @ one gets 

a = a Q t^TV 7 (6) 

or, in other terms, 

h = WTWf (7) 

Another useful equation obtained by differentiation of the Eq. © and substitution into it Eq. reads: 

h = -\{ P + p). (8) 

The pressure can be written down as 

P=-\h-h 2 . (9) 

It is well known that for a given cosmological evolution h(t) satisfying some simple conditions one can find a 
minimally coupled scalar field cosmological model with Lagrangian 

L=¥--V(<f>), (10) 

which contains this evolution as a particular solution, provided suitable initial conditions are chosen (see, e.g. |l0lll9| 
and references therein). 

The energy density and pressure of the scalar field 1)10(1 are 



p=\-v(4>). (12) 
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From Eqs. (HJ, JT^J, © and (SJ one has 



2 = (p + p) = -|ft, (13) 



^) = ^-p) = ~ + /* 2 - (14) 

Equation (fT4"|) gives the potential V(t) as a function oft. Finding </>(i) by integrating Eq. lfT3)l. inverting this function 
and substituting a function i(0) into Eq. JT2J, one finds the potential V as a function of the scalar field 4>. Sometimes 
this can be done explicitly. For example, the evolution Q could be reproduced in the extensively studied 0, |2(j 
cosmological model with an exponential potential 



n*)~«p(-*£p*Y as) 



Apparently, the super-accelerated evolution Q with fc < — 1 cannot be reproduced by means of the minimally 
coupled scalar field model 1)10(1. but can be easily obtained in the framework of the model with phantom scalar field 



£ = -y-W (16) 



Our starting point is the remark made in |21|, l22t |23j , attracting the attention to some interesting features of the 
potential i|14|) represented as a function of time t. The point is that the volume function 

m = a 3 (t) (17) 

satisfies a simple second-order differential equation 

$ = 9V(t)il>, (18) 

which can be easily obtained from Eq. I|14fl . As has been already noticed above, the potential V(t) is chosen in such 
a way to provide a given cosmological evolution. For example, for the evolution the form of the potential is 

2(1 -fc) 

V ^ = 9(TTW- ^ 

For simple forms of the potential one can find the general solution of Eq. (|18fl . which contains together with the 
solution used for the construction of this potential also another independent solution. For the potential ijT§|) the 
general solution looks like 

m = ^it ai + ^t a \ 

2 

ai = TTI' 

« 2 = rri (20) 

where rj}\ and -02 are nonnegative constants. (In papers [2^, |2^| the solution l|2(JI) was written down for the case of 
dust k — 0.) Now, knowing ip(t) one can find h(t) = A and h. Substituting the obtained h(t) into Eq. I|14l) one 
can try to find the dependence <j>(t). Inversion of this dependence t(4>) and its substitution into the expression for the 
potential as a function of time l|19|) would have given the potential V(4>), whose form is different from the exponential 
one (fl3|l when \jj 2 ^ 0. 

One can derive the equation for the functions t(4>) compatible with the potential V(t) from Eq. I|15fl. First of all, 
notice that the field cf> should satisfy the Klein-Gordon equation 

+ 3/10+^ = 0. (21) 
dep 



Dividing Eq. (|21(l by cb one has 



3/i+i^ = 0. (22) 
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It is convenient to rewrite the second time derivative of the scalar field as 
where we have used the relation 

f = i, (24) 

with the "prime" denoting the differentiation with respect to the scalar field <j>. Now, we would like to take the time 
derivative of Eq. The first term gives 



1m = L(1 

dt\d>) f\f 



(25) 



Then, using the relations l|14|l and (|24[1 one gets 



Then we make the following transformation: 



l^r=t"V. (27) 



The time derivative of the expression in the left hand side of Eq. (|27|l is 

^(t 2 V)'. (28) 



dt 



Combining expressions l|25 |l .l|26 |l and (|28|l . using the explicit form of the potential Ijl9|) and multiplying the resulting 
equation by t' one arrives to the following equation: 

1\" 9 4(1 -k) ( t" ' 



w 9(i + kyyw) =°- (29) 

Particular solutions of this equation being substituted into Eq. H19|l generate different potentials V{<f>) corresponding 
to different cosmological evolutions (|2U|I . It is easy to check that the solution 



t = exp - , (30) 

corresponding to the potential l|15|) and to the evolution with ip 2 = 0, satisfies Eq. (|29|l . Instead of looking for other 
exact solutions (which is not an easy task) we concentrate now on the qualitative study of Eq. I|29|l or its equivalents. 

First of all, let us notice that the solution l|20(l tacitly assumes that the time parameter runs from to +oo. When 
t — ► +oo one can neglect the influence of the second term in the solution 12Ufl . Thus, the potential V(4>) in this limit 
has the asymptotic form 115fl . It is possible to calculate also the subleading term. The form of the Hubble variable 
corresponding to the evolution (|20|l is 

, ai^ii" 1 - 1 +a 2 ^ 2 P 2 - 1 

n — ; ; , (-jl 

3(-0ii Ql +V>2^ 2 ) 



while its time derivative is 



• _ a^jt 2 ^ + a 2 y.|t 2Q » + (1 - (ai - a 2 ) 2 )^ 2 t a ^ 
1 3i 2 (V>ii Ql + V^" 2 ) 2 ' 



At large values of t, h written down up to the first non-leading term reads: 

2 ( 2(fc-3)V 
3(l + /c)i 2 ^ (l + fc)V 



2 / 2(^-3)^2 *-»\ 
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Substituting the expression (|33[1 into Eq. I|14|) one finds the subleading correction to the scalar field <f> as a function 
of t. Inverting this relation and inserting the result into Eq. I|19|) we obtain the potential V((f>) including the first 
nonleading term: 



i- f ' -3vT + ffi v / , V> 2 fc-3 3(fc-3)^ 
V{(p) ~ exp ( 1(1 + 3- 7= ; exp 



VTTfc 



2\/T 



(34) 



The above formula can be elucidated as follows: if one consider the cosmological evolution described by Eq. i|2(J|) 
with non-zero value of the parameter ip2, one can state that there exists a potential describing this evolution. In spite 
of the absence of an explicit expression for this potential we know that its asymptotic behaviour is described by the 
formula 

Now, let us consider the evolution of a cosmological model, which is characterized by Eq. (|31|l and the potential 
whose asymptotic form is given by Eq. H34|) . It is easy to show that at some moment to, such as < <o < 00 , the 
time derivative of the Hubble variable vanished. Indeed, a simple analysis of the formula (|32|) gives this value 



(3 - k)y/2(l-k) + 4(1 - k) V>2 



2(1 + A) 



(35) 



Obviously, to exists for any value of ip2 0. The evolution of the model given by Eq. I|31(l is consistently defined at 
t > to in the standard scalar model (|10[1 with a potential whose asymptotic form at large values of t was presented 
above H34fl . (The consistency is based on the fact that at t > to the time derivative of the Hubble constant is always 
negative). It is possible also to find the approximate expression for the potential V(cf>) for t — ► to+- it is well known 
that the Klein-Gordon and Friedmann equations allow a shift of the scalar field by some constant. It is convenient to 
fix the scalar field as 4> — > when t — > t +- The time derivative of the Hubble variable at t — > t + behaves as 



where 



H 



h = -H(t-t ), 



y/8(i-k)(3 - k)H'ii'2t 14 



(i + fe) 3 +MT 1 



(36) 



(37) 



Substituting the expression H3HJI into Eq. I|14|l one can integrate the latter, getting 




(t-tof 2 , 



whose inversion gives 



Thus, the potential results 



t = t + 




V{cj>) 



2(1 - k) 




(38) 



(39) 



(40) 



and the potential V(<j>) are regular when t — > fo- 



rt is rather straightforward to verify that the expressions for 
while the derivative dV J^ and <j> are singular. 

Now, before discussing the subtle question of the behavior of the cosmological model under consideration at t = to 
and of the possibility of crossing this point, we shall study briefly another branch of the evolution l|31l) : namely, 
that where t < to- For small values of t to reproduce the evolution (|2C)|> one should use the phantom Lagrangian 
()16J) . Proceeding in the similar way one can find the corresponding potential of the phantom field in the subleading 
approximation: 



t/y j\ 2(1 -k) / 2(1 + k) 



1 



2(1 -k 2 



(41) 
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One can easily describe what is going on in the model characterized by a negative kinetic term and the potential 
whose asymptotic form is given by (|41|l in the vicinity of the point t — > t — . 

Now, we are in a position to scrutinize the possibility of matching these two branches of the evolution 1)20(1 at t < to 
and t > to- Are they really incompatible? Let us consider a scalar field model with a negative kinetic term and a 
potential whose asymptotic form is that of Eq. (|41|l . Further, let us suppose that at some moment t < to one has 
initial conditions on the values of h, <fi and <p which provide the evolution (|2(J|) where the values of the parameters are 
consistent with those of the potential V(4>). Then approaching the time moment t — > to— one arrives to the regime 
when h, <fi and <f> tend to vanish. Nevertheless, all the geometric characteristics of the spacetime remain well defined 
(due to homogeneity and isotropy of the Friedmann cosmology all the curvature invariants are expressible through 
the cosmological factor a and its time derivatives, which are obviously finite). In contradistinction, the second time 
derivative of the scalar field at t = to and dV J^ diverge, but this divergence is an integrable one. This offers us an 
opportunity (and, perhaps, necessity) of a continuation of the spacetime geometry and field configurations beyound 
this "divide line" . Clearly, such a smooth dynamical continuation which respects Einstein equations, entails the 
change of the sign of the kinetic term and the transition to the regime for t > t described above. 



III. DIFFERENT REGIMES OF CROSSING OF PHANTOM DIVIDE LINE 



The cosmological evolution presented in the preceding section describes the following scenario: the universe begins 
its evolution from the cosmological singularity of the "anti-Big Rip" type and its squeezing is driven by a scalar field 
with a negative kinetic term. Then at the moment (easily calculated from Eq. 

/(l-fc)^A^ , , 

the contraction of the universe is replaced by an expansion. At the moment t = t > t\ (see Eq. I|35[l the kinetic term 
of the scalar field changes sign. Then, with the time growing the universe undergoes an infinite power-law expansion. 
The graphic of the time dependence of h(t) is presented in Fig. 1. Changing the sign of time parameter in all the 




FIG. 1: h(t) dependence in our toy model 



equations the range — oo < t < one can consider the cosmological evolution running from t = — oo to t = 0. It 
begins from the non-singular contraction and ends in a "standard" Big Rip singularity. In this case at the moment 
t = —to the scalar field undergoes the "phantomization" , while the transition at the moment t = to, considered in the 
preceeding section can be called "dephantomization" . 

The toy model discussed above prrovides an illustration of these phenomena of phantomization and dephantomiza- 
tion, but from observational point of view, it would be more interesting to get an evolution beginning from the 
Big Bang and ending in the Big Rip singularity, after undergoing a phantomization transition. Instead of trying 
to construct some potential and cosmological evolution describing such a process we shall limit ourselves to giving 
a graphical presentation of the h(t) -dependence which could be responsible for such a scenario (see Fig. 2). This 
picture represents the function h(t) — t ( t ^-t) ' Here the moment t = tR is the moment of Big Rip, while the the 
moment of phantomization is to = tu/2. In principle, the function h(t) determines the model completely, including 
the form of the potential, even if when it is not written explicitly. 

As has been already mentioned in Introduction, some observations favor double crossing of the phantom divide line 
. The corresponding graphic h(t) is presented in Fig. 3. Also to this time dependence of h(t) one can associate some 
potential and the evolution, in course of which the scalar field undergoes two smooth transitions: phantomization at 
the point of mimimum of the curve h(t) and the subsequent dephantomization at the point of its maximum. 



FIG. 2: h(t) dependence in the model, describing the Big Rip. 



h[t] 



FIG. 3: h(t) dependence in the model, describing a double crossing of the phantom divide line. 



It is suggestive to remark that the diagram Fig. 3 reminds topologically the well-known van der Waals isoterm curve. 
The part of the van der Waals curve situated between the maximum and minimum, describing a metastable state (of 
supercool vapor or superheated liquid), is analogous to the phantom pha se in our Fig. 3. Models with double crossing 
of the phantom divide line were considered also in p2^ . The idea in |22j consisted in cutting off the phantom piece of 
the evolution with the subsequent sewing of remaining branches of the trajectory removes the problems connected with 
the crossing of the phantom divide line. In other terms, the trajectories, considered in p2^ correspond graphically to 
h(t) which instead of one maximum and one minimum (sec Fig. 3) have one inflection point. Continuing the analogy 
with the van der Waals isoterms, one sees that this curves corresponds to the critical van der Waals curve. Surely, 
such models can be studied, but their existence does not exclude the existence of models described by Fig. 3 and in 
a frame of a concrete model with fixed Lagrangian and initial conditions these two types of trajectories cannot be 
transformed one into other. 



IV. DISCUSSION 



Above we have discussed the properties of some simple cosmological models based on minimally coupled scalar 
fields, which could, in our opinion, describe the crossing of the phantom divide line, changing the sign of the kinetic 
term in the corresponding Lagrangian. As a starting point of our consideration we have used the observation |22ll23| 
that having a simple form of the scalar field potential , written as a function of time V(t), one can find the general 
solution of the Einstein (Friedmann) equations for the Hubble variable h(t) (see Eq. (J2DJ. Different particular solutions 
describing the evolution of the Hubble variable correspond to different potentials as functions of the scalar field V(cf>). 
All these solutions (excluding the obvious, when tp2 — 0) contain the time moment when the time derivative of the 
Hubble variable vanishes (see Eq. ifSBfr) and the only consistent way of passing of this point respecting the continuity 
of the Einstein equations is the change of the sign of the kinetic term for the scalar field. This change can be called 
(dc)-phantomization of the field. 

All above are of technical nature and are confirmed by a direct analysis of the corresponding differential equations. 
As far as physical interpretation is concerned one has two alternatives. One can say that the sign of the kinetic term 
should be fixed a priori (see, e.g. [lj|)) and, hence the models, which imply the change of this sign should be discarded 
as non-physical ones. Here, we attempted to treat another alternative: i.e. to take this theory seriously because from 
mathematical point of view it looks consistent. In choosing this alternative we were guided by a belief that Einstein 
equations are more fundamental then the concrete form of the action for other fields. The idea about the dominant 
role of Einstein equations goes back to the classical works by Einstein, Infeld and Hoffmann pi. l25| , where it was 
shown that the motions of mass points (geodesies law) are determined by the Einstein equations for the gravitational 
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field. Later this approach was confirmed in works by Lanczos [2(j, Fock [23 and others. In this respect the general 
relativity is quite different from electrodynamics, where the Lorentz's force law is quite independent of the Maxwell 
field equations (for a discussion see e.g. Hi])- The derivability of the equation of motion for particles is based on the 
fact that the Einstein equations are non-linear and also they are subject to additional identities, namely the Bianchi 
identities which reduce the number of independent equations. These Bianchi identities imply conservation laws for 
matter, present in the model under consideration. 

In fact, the argumentation, based on the Bianchi identities could be used also for getting information about the 
possible field configurations for non-gravitational fields. Indeed, the Einstein equations for fields, minimally coupled 
to gravity, have the form 

K - = kT£. (43) 

The Bianchi identities 

v M (m - \kr) = (44) 

imply some kind of the energy-momentum conservation law: 

= 0. (45) 

For the case of the minimally coupled spatially homogeneous time-dependent scalar field with the standard sign of 
the kinetic term Eq. Q45JI is reduced to 

+ 3/i0+^)0 = O, (46) 



3/i0+^)0 = O. (47) 



which is equivalent to the Klein-Gordon equation when ^ 0. For the phantom scalar field one has similarly 



The comparison between Eqs. <|46[) and (|47|l points to the opportunity of the change of type of the Klein-Gordon 
equation at the moment when = 0. In other words, while the Einstein equations require the change of the type of 
the scalar field Lagrangian provided some form of the scalar field potential and the cosmological evolution are chosen, 
the Bianchi identities show why this transformation is possible. 

Before concluding, let us notice that we have shown that simple scalar field models driving the universe through 
the phantom divide line do exist. However, the natural question arises: given the potential as a function of the scalar 
field (not of the time) how general are the initial conditions providing the reaching of the point = 0. In other words, 
do we need the fine tuning (cf. [l8[ l ?. The answer to this question is not simple technically and, perhaps, is worth 
further investigations. Another interesting question is concerned with an opportunity of using some kind of modified 
local Lagrangian formalism for the description of smooth dynamical transition from one form kinetic term to another 
one. 
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